Curvature effects on the surface thickness and tension at the free interface of 4 He 

systems 
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The thickness W and the surface energy o~a at the free in- 
terface of superfiuid 4 He are studied. Results of calculations 
carried out by using density functionals for cylindrical and 
spherical systems are presented in a unified way, including 
a comparison with the behavior of planar slabs. It is found 
that for large species W is independent of the geometry. The 
obtained values of W are compared with prior theoretical re- 
sults and experimental data. Experimental data favor results 
evaluated by adopting finite range approaches. The behavior 
of a a and Wo a exhibit overshoots similar to that found pre- 
viously for the central density, the trend of these observables 
towards their asymptotic values is examined. 

PACS numbers: 61.20.-p,68.10.-m,68.10.Cr 



I. INTRODUCTION 

The understanding of the density profiles in the sur- 
face region of a quantum fluid like 4 He has long been 
considered a very important basic problem. 1-3 At the 
liquid-vapor interface, the density profile of 4 He changes 
continuously from liquid density pi to vapor density p v 
over a distance of some angstroms. In the case of a 
liquid-vacuum interface at T = K, pi falls monoton- 
ically to p v = 0. The width W of a surface is defined 
as the distance in which the density decreases from 0.9p£ 
to 0.1p£. A glance at recent literature reveals several 
works addressing the question of the thickness of the 
free interface, 4 " 7 indicating the continuous interest in this 
area of theoretical and experimental research. The sur- 
face tension at the free interface has been also investi- 
gated for a long time. 1 ' 5,8-12 A list of results for the sur- 
face thickness and tension determined up to the middle 
of 1987 is given by Osborne in Table 1 of Ref. 2. 

In the systematic study of free planar 4 He films at 
T = K we have, among other issues, discussed features 
of the surface thickness. 5 There, our own results for W 
evaluated by using several density functional (DP) ap- 
proaches are compared with values obtained from Monte 
Carlo simulations by Valles and Schmidt 12 and experi- 
mental data of Lurio et al. 3 From Fig. 4 in Ref. 5 one 
can realize that the size of the experimental error bar 
was to large for disregarding any of the applied DF ap- 
proaches. After that paper, new theoretical and experi- 



mental results for this quantity have appeared. Evalua- 
tions of W for free planar slabs and droplets by utilizing 
a variational Monte Carlo with shadow wave functions 
with a glue term (glue-SWF) were published by Galli 
and Reatto. 6 On the other hand, a novel measurement of 
W by using x-ray reflectivity was reported by Penanen et 
al. 7 These new data supersede that previously obtained 
at the same laboratory. 3 To complete the survey of inves- 
tigations about the surface thickness of superfiuid helium 
droplets, one should mention the calculation of Stringari 
and Treiner 13 and the comprehensive experimental and 
theoretical study of Harms, Toennies, and Dalfovo. 4 

The foregone summary indicates that there is an im- 
portant piece of information about the free surface thick- 
ness of planar and spherical superfiuid 4 He systems. Al- 
though in recent years there is a renewal of interest for 
examining cylindrical species, 14-20 hitherto, there is no 
study undertaking the problem of exploring features of 
W in the case of such geometry. 

On the other hand, we have explored the evolution of 
the surface energy at the liquid-vacuum interface of pla- 
nar slabs as a function of their size in Ref. 5. However, 
as far as we know, there is no work devoted to study 
this property in the case of systems with curved geome- 
tries. Therefore, such an analysis becomes an interesting 
problem by its own right. 

In view of the situation described in the previous para- 
graphs, the aim of this work is to study systematically the 
interface thickness and the surface tension of liquid 4 He 
with cylindrical and spherical shapes making a connec- 
tion to the case of planar slabs. Although the spherical 
systems are energetically favored against free cylinders 
and slabs an analysis of the overall picture presents in- 
structive features. The width at the liquid-vacuum in- 
terface of free systems is compared with experimental 
data and with theoretical results obtained for stable pla- 
nar films of 4 He adsorbed onto the lightest alkali metals. 
The theoretical tools are outlined in Sec. II. In Sec. 
Ill we search for the size at which the systems reach an 
asymptotic global behavior. The discussion of the pat- 
tern exhibited by the width and tension at the surface 
is done in Sec. IV, where the results are presented in 
a unified way allowing a direct comparison of data ob- 
tained for different geometries. Section V is devoted to 
a summary. 



1 



II. THEORETICAL FRAMEWORK 



A. Zero-range density functional 



The calculations performed in the present work were 
carried out by using density functional (DF) approaches, 
which have proven to be successful tools for treating this 
kind of quantum many-body problems. In such a theory 
the ground-state energy, E gs , of an interacting TV-body 
system of 4 He atoms, confined by an adsorbate-substrate 
potential C/ su b(r), may be written as 

E gs = J dvp{v)H[p,Vp) + J dr p(r) U snh (r) 

= -|^ / ^VpWV 2 v / pW+ J drp(r) e sc (r) 
+ fdrp(r)U suh (r) , (2.1) 



where p(r) is the one-body density. The first term on 
the right-hand side is the quantum kinetic energy of the 
helium particles of mass m. The second term represents 
the interaction between the particles of the system, where 
e sc (r) is the self-correlation energy per particle depending 
on the DF approach. The last term is the interaction with 
the external field. 

The density profile p(r) is determined from the Eulcr- 
Lagrange (EL) equation derived from the condition 



Sn _ 5{E gs [p,Vp}- pN} 



6p(r) 



Sp(r) 



= . 



(2.2) 



Here p is the chemical potential, N the number of parti- 
cles 



N 



j dvp{r) , 



(2.3) 



and £1 the grand thermodynamic potential. The variation 
of (2.2) leads to a Hartree like equation for the square 
root of the one-body density 



2m 



V 2 + V H {v) + U suh {v) 



VpW^VpW, (2-4) 



which also determines p. Here Vjj(r) is a Hartree mean- 
field potential given by the first functional derivative of 
the total correlation energy E sc [p] 



Mr) 



6E sc [p, Vp] 



5p(r) 
d 



- V 



Sp(r) 
d 



dp(r) dVp(r) 



J dTlp{vl) e sc (r/) 
J dvi p(ri 



)e sc (r/). (2.5) 



Two different DF were used, namely, the Skyrme type 
"zero-range" version suggested in Refs. 13 and 21, and 
the nonlocal density functional (NLDF) proposed in Ref. 
22. For each DF an expression for Vn(r) should be de- 
rived. 



The simplest DF successfully employed to interpret 
properties of 4 He systems is a zero-range correlation pro- 
posed by Stringari and Treiner. 13 ' 21 It has been inspired 
in Skyrme type functionals extensively used to describe 
properties of atomic nuclei. 23 ' 24 The explicit form of this 
correlation energy per particle is 



4 k5 » = 



P(r) + yP 74+1 (r) + d 4 



1 



p(r) 



Vp(r) | 2 . 
(2.6) 



The phenomenological parameters 64, c 4 , 74, and d 4 have 
been fixed in Ref. 13 so as to reproduce the known ob- 
servables of the bulk liquid at equilibrium. The data of 
these saturation quantities (where P = 0), i.e., the equi- 
librium density po, the minimum energy per particle es, 
the compressibility /C, and the surface tension of a 
semi-infinite 4 He system are listed in Table I. Experi- 
mental values are correctly reproduced by the set 



b 4 = -8.888 10 x 10 2 KA 3 , 
c 4 = 1.045 54 x 10 7 KA 3 <^ +1 ) 
74 = 2.8 , 

di = 2.383 x 10 3 KA 5 . 



(2.7) 



B. Orsay-Paris nonlocal density functional 

Currently, in the literature one may find a few NLDF 
approaches. The Orsay-Paris (OP) version (OP-NLDF) 
developed by Dupont-Roc et al. 22 treats correctly the 
long-range part of the helium-helium interaction and pro- 
vides a reasonable description of correlations. The most 
elaborated version has been formulated by the Orsay- 
Trcnto (OT) collaboration. 25 However, the OP-NLDF is 
sufficiently good to reproduce properties of non-layered 
samples like free or weakly confined systems. This func- 
tional reads 



.op 



(r) 



-V 



c2r/p(r/)V? JF (| r - r/ 



c 4 



[p(r) 



174 + 1 



(2.8) 



OP/ 



r — r/ 



In this case the two-body interaction, V l 
was taken as the 4 He- 4 He Lennard- Jones (LJ) potential 
screened in a simple way at distances shorter than a char- 
acteristic distance hop 



Vr(r)={ 



4e LJ 



\ r ) V r J 



if r > hop 



if r < h p 



(2.9) 



with the standard de Boer and Michels parameters, 26 
namely, well depth 6l.t = 10.22 K and hard core radius 



2 



(Tlj = 2.556 A. In order to recover the correct results 
for bulk liquid, the screening distance hop was adjusted 
so that the integral of V t (r) over the whole three- 
dimensional space be equal to the value of 64 given by 
Eq. (2.7) 



b 4 = J drV^ir) 



32tt 
"21 



■^LJ £ LJ 



8 ( <TLJ 



3 \hop 



3n 



(2.10) 



This procedure led to /lop = 2.376 728 A. 

The p(r) is the "coarse-grained density" defined as the 
straight average of p(r) over a sphere centered at r and 
with a radius equal to the screening distance hop 



(2.11) 



p(r) = / dvi p(vi) W(\ r - r/ |) 



where W(| r — r/ |) is taken as the normalized step func- 
tion 



W (|r-r/|) = i -^e(^ p-|r-r/|) 



■-op 
3 



if I r - r/ I < h p , 



(2.12) 



if I r - r/ I > h p ■ 



a rather simple compact form similar to Eq. (2.16) writ- 
ten in Ref. 27 for planar films instead of adopting the 
expansion in terms of Legendre polynomials proposed in 
Ref. 28. The latter procedure is more appropriate for 
studying excitations of a given system. We shall mainly 
report results for curved geometries obtained by setting 
U suh (r) ee 0. 

In the case of cylindrical symmetry Eq. (2.13) is solved 
for a fixed number of particles per unit length L 



n\ = N/ L = 2 7r / r dr p(r) 
Jo 



(2.15) 



while for helium spheres the constrain is directly the 
number of particles 



N 



: 7T / 

Jo 



dr p(r) 



(2.16) 



Typical density profiles are displayed in Fig. 2 of Ref. 18 
for cylinders and in Fig. 1 of Ref. 13 for spheres. 

We shall first test whether the size of the considered 
systems of liquid He is large enough to expect asymp- 
totic behaviors. Next, we shall focus our attention on the 
surface energy and the width of the interface. 



III. ASYMPTOTIC GLOBAL BEHAVIOR OF 
THE SYSTEMS 



C. Hartree like equation 

In the case of curved geometries the Eq. (2.4) takes the 
form 

h 2 ( d 2 D - 1 d \ r-r-: 



+ 



V H (r) + U suh (r) 



y/W) = V JpVj , (2-13) 



with D = 2 for cylindrical systems and D = 3 for spheri- 
cal ones. The equation for planar systems is obtained by 
setting D = 1 and assuming that r represents the coor- 
dinate perpendicular to the plane of symmetry (usually 
denoted z). The Hartree potential derived by applying 
Eq. (2.5) to the Skyrmc-DF reads 

V^(r)=b 4 p(r) + 2±±l CiP ^+\r) 

„ , / d 2 D- 1 d\ . . , n A , 

The expression for Vu(r) derived in the OP-NLDF ap- 
proach for cylindrical systems is given in the Appendix 
of Ref. 18. The corresponding one for spherical systems 
is provided in the Appendix A of the present work. It is 
worthwhile to notice that we expressed this quantity in 



In the case of spherical systems, it has become custom- 
ary to define for each moment < r k > an equivalent uni- 
form radius Rk- This quantity is equal to the radius of a 
uniformly occupied sphere (of density p u ) with the same 
momentum of order k given by the true density distribu- 
tion. It is possible to extend this idea to the cylindrical 
geometry. So, in general, an equivalent uniform radius 
Rk associated to a momentum r k given by 



2 D-i nL 3-D jRk dr r k+D-i pu 



<<>D = 



n L*~ D f* k dr r D ^ p 



J Rk dr r k + D ^ p v 



_ Jo 



D 



f" k drr D -^p u k + D 



R\ , (3.1) 



is determined by matching the result of this integral with 
the k- moment of the true radial distribution of the system 



<r k (N) > D - 



2 D - 1 7rL 3 - D 
N 



f 

Jo 



dr r k 



k+D-l 



p{r) . (3.2) 



This condition leads to the generalized Ford- Wills mo- 
ments (see Eq. (1) in Ref. 29) 



Rk 



D 



D 



< r k (N) > D 



l/fc 



(3.3) 



3 



A. Cylindrical systems 

One expects that very thick cylinders should tend to 
exhibit features of bulk liquid at saturation conditions. 
In the asymptotic limit of very large cylinders one expects 
that Eq. (2.15) will reduce to 



N/L = 2 7r / r dr po = n Rq po , 



(3.4) 



where po is the saturation density of bulk 4 He quoted in 
Table I and Rq is the radius of the uniformly filled cylin- 
der. Since in the limit of very large cylinders the equiva- 
lent radius Rk tends toward Ro defined by Eq. (3.4), one 
gets the following asymptotic law 

i/fc 



lim < 

ra A ->oo 



[< r k (n x ) >]V* 



n 



1/2 



1 

7T/90 



k + 2 



(3.5) 



Figure 1 shows how the normalized moments < r > and 
< r 2 > attain their asymptotic values as a function of the 



expansion parameter v cy \ = n 



-1/2 



One may conclude 



that the results given by Eq. (3.5) are safely reached for 
n A ~ 70 A" 1 . 



B. Spherical systems 

From a similar analysis to that performed for cylinders, 
one can state that for very large spheres Eq. (2.16) will 
become 



N 



Jo 



2 j 4t „•! 

r z dr p = — i?g po , 



(3.6) 



where Rq is now the radius of the uniformly filled sphere. 
Furthermore, since in this asymptotic limit Rk tends to- 
ward R defined by Eq. (3.6), the following law is ob- 
tained 

1/3 / o \ 1/fc 

lim 



N- 



[< r k (N) >} 1 / k \ _ / 3 
N 1 / 3 J ~~ V 47r /°o 



fc + 3 



(3.7) 

Figure 2 shows how the normalized moments < r > and 
< r 2 > attain their asymptotic values as a function of 
^s P h = iV -1 / 3 . The results given by Eq. (3.7) are reached 
at N ~ 2000. 

In summary, from values of the normalized moments 
displayed in Figs. 1 and 2 it is clear that the largest 
systems examined in the present work have, in practice, 
reached the asymptotic global behavior. 



IV. ANALYSIS OF RESULTS 

A. Surface tension 

The evolution of the central density, p c , for cylinders 
with increasing n\ has been explored in Ref. 18. It was 



found that after a certain value of n\ the central density 
becomes bigger than the saturation density of infinite 
helium matter, p c > po, giving rise to a squeezing effect. 
This phenomenon has been previously found in calcula- 
tions carried out for atomic nuclei (see Figs, la and lb 
in Ref. 30) as well as for spherical clusters of 4 He (see 
Fig. 3 in Ref. 13) and it is known as the "leptodermous" 
behavior. It was demonstrated 18 that for cylinders the 
squeezing effect, i.e. the difference p c — po, vanishes for 
n\ — > oo according to the law 



e cy i = 

which is similar to 

Pc - Po 




e sph 



PO 



,3PoJ 



K 



(4.1) 



(4.2) 



corresponding to 4 He spheres (see the asymptotic limit 
of Eq. 12 in Ref. 13). 

In the present work we shall describe the evolution of 
the surface tension when the size of the He systems is 
increased. From thermodynamic considerations one gets 
that 



dE gs = -PdV + o A dA + pdN 



(4.3) 



where V is the volume of the system, A the surface area 
enclosing V, and a a is the surface tension. For both 
curved geometries the grand free energy takes the form 



O = E gs - p N = -P V + a A A , 



(4.4) 



while P = for planar slabs. In the following lines we 
shall examine the different cases. 



1. Cylindrical systems 

In Ref. 18 it has been shown that for cylinders the re- 
lationship between the grand free energy per unit length 
and the surface tension becomes 



Q E gs — pN 



(e-p)nx =ttR s <7 A ■ (4.5) 



L L 

Here R s is the sharp mean radius defined by 

N/L = 2 7T / r dr p c — tt R 2 p c , 
Jo 

and it is related to Rq by 



Rs — Ro 



1 



1 + Ccyl 



(4.6) 



(4.7) 



Furthermore, it has been demonstrated in Ref. 18 that 
for large cylindrical systems (in the sense of big R s ) the 
total energy per particle may be expressed as 



4 



e = 



N 




2p a lC 



(4- 



where Coo and ee are, respectively, the asymptotic surface 
tension and the residual energy per unit length defined 
in Ref. 18. A formula for the surface tension a a may be 
derived by taking into account Eq. (5.7) in Ref. 18 



D 2 de 

ttR s a A = -n x - — 



(4.9) 



and the present relationship (4.5). After keeping all 
ip to 2-nd order in i 

2ft _ 2{e-p,)N 



terms up to 2-nd order in v cy \ = n x one gets 



a a = 




de 



it R s dn\ 



3 



£t. 7 



4 p JC 



(4.10) 



This expression indicates that the asymptotic result 
oa(it>\ — > oo ) = doc, which may be identified with 
the experimental 8-10 values of the surface tension er cxp 
quoted in Table I, should be attained linearly in the pa- 

— 1/2 

rameter v cy \ — n x with a slope proportional to q 



Po -1/2 

o A ^ o-oc + \ —e e n x 

TT 



(4.11) 



The values of Ei obtained from the fits of energies per 
particle to Eq. (4.8) are listed in Table I. Since these re- 
sults are positive one expects an overshoot (i.e., a region 
where (j A {n\) > (Too) for large cylinders. 



with R defined by Eq. (3.6). Then by using (4.4) one 
gets 



ft = -PV + a A A = - 



2 a a 

Rs 



(4.16) 



The total energy per particle can be expanded in the 
following way [cf. Eq. (13) in Ref. 13] 



N 



e B + (J c 



36 it 



Pi 



1/3 



+ 



2^(^\ 



2/3 



N - 1/3 

N~ 2/3 + --- . (4.17) 



Hence, for spherical systems one gets 



OA 



3ft _ 3(e-/x)AT 



3 N 2 de 



^4sph 

o"oo + a c 



AttR 2 



AttRI dN 



±Po\ N -i/ 3 
9tt 



(4.18) 



The values of a c determined by fitting energies per par- 
ticle of drops with N > 300 to Eq. (4.17) are quoted in 
Table I. Since these results are positive, overshoots of 
a A should be expected. For the sake of comparison prior 
results reported in Refs. 13 and 21 obtained by using 
the Skyrmc-DF and considering systems with N < 728 
are also included in Table I. All these data are in good 
agreement. 



3. Comparison of results for regular geometries 



2. Spherical systems 

In the case of spherical systems, upon starting from 
Eqs. (4.3) and (4.4) one arrives at 

dft = d(E gs - pN) = -PdV + cr A dA- Ndp . (4.12) 

At equilibrium, for a fixed /x, the virtual work by chang- 
ing the radius of the sphere should vanish 

-PdV + a A dA = -P4:TrR 2 dR s +<T A 8irR s dR s = . 



This equation leads to 



2 OA 



where the sharp mean radius is 
Rs = Ro 



1 + e sp h 



1/3 



(4.13) 



(4.14) 



(4.15) 



Let us now summarize the results for the surface ten- 
sion for the three regular geometries. Instead of show- 
ing the surface tension separately for each geometry as 
a function of the appropriate expansion parameter v 1 we 
find it more interesting to present all the data in a uni- 
fied picture. In order to achieve this goal it is useful to 
consider that for all three symmetries (planar, cylindrical 
and spherical) the final trend of the energy per particle 
toward the asymptotic value, i.e., the bulk es, is mainly 
determined by the linear term in the corresponding ex- 
pansion parameter. Therefore, it becomes reasonable to 
adopt the energy difference e — es as an appropriate com- 
mon abscissa. In order to facilitate the forthcoming dis- 
cussion we recall that for a symmetric planar slab (see, 
e.g., Ref. 5) the energy per particle may be expressed as 
an expansion in terms of inverse coverage n^ 1 — A s i a b/A 



E 

e = -^p = e B + 2 (Too n- 1 + 7c n -3 + 
and the total surface tension becomes 



(4.19) 



•5 



(tot) 



E gs - pN _ (e- p)N 



o de 



= - n r 



L 2 

2 (Joo + 3 7 C n~ 2 . 



(4.20) 



Values for the surface tension obtained from calcula- 
tions carried out by using the Skyrme-DF are displayed 
in Fig. 3, while Fig. 4 shows the behavior of a a evaluated 
by using the finite range OP-NLDF. The main features 
exhibited by both these drawings are similar. In partic- 
ular, for curved geometries the surface tension tends to- 
ward (Too showing the predicted overshoot, which is sim- 
ilar to that previously found for p c (see Fig. 6 in Ref. 
18 and Fig. 3 in Ref. 13). Furthermore, in both figures 
the final trends toward match the linear asymptotic 
expressions: 

a) for cylindrical systems 



a A — o"o 



Po £i 

2 7T (Toe 



b) for spherical systems 



(e - e B ) ; 



^(^Y /3 (e-e B) . (4.22) 

3 7T (Too V b / 



(4.21) 



(TA 



On the other hand, the data for planar slabs attain the 
asymptotic value in a smooth way with zero slope, that 
is without any linear term. This behavior may be un- 
derstood on the basis of the fact that upon starting from 
Eq. (4.20) one gets 
c) for planar systems 



(tot) 



OA 



— COO + 



3 7c 
8^ 



(e - e B f 



(4.23) 



Here a a is the surface tension at each free interface. The 
values obtained with the OP-NLDF approach match this 
law with 7 C = - 1 .44 x 10~ 3 K /A 6 listed in Table 2 of Ref. 
5, while the trend of results provided by the Skyrme-DF 
is even flatter suggesting a (e — es)" law with n > 3. The 
latter behavior is due to the fact that in such a mean- 
field approximation j c ~ because the Lennard-Jones 
potential is not included in the theory (cf. the Table 2 
and the discussion in Ref. 5). 



4- Relation between a c and ee 

We shall now search for a relation between a c and 
e e - According to the Droplet Model (DM) 30 ' 31 the to- 
tal ground-state energy of a large system may be written 
as a sum of volume, surface, and curvature terms 



i?g S — E v + E s + E c . 



(4.24) 



The surface and curvature energies for cylindrical sys- 
tems may be written as (see Eq. (4.8) in Ref. 18) 



poo 

E s (cy\) = 2nR L {p(r)H[p,V p] - p Q e a }dr 
Jo 

= 2^i? (cyl)ia co , (4.25) 

£ c (cyl) =2ttL / {p(r)H[p,\7p}-p e }(r-R )dr 
Jo 

= s e L, (4.26) 
while for spheres according to Eq. (2.11) in Ref. 31 holds 

/>oo 

£ s (sph) =4 7ri?2 / {p( r )H[p,Vp} - Po e }dr 
Jo 

= AnRl(s V h)a oc =a s N 2 / 3 , (4.27) 

E c (sph) =8irR / {p(r)H[p,Vp} - p Q e Q }(r - R )dr 
Jo 

= a c N 1 / 3 . (4.28) 

In this approach the integrals in Eqs. (4.25) and (4.27) 
are associated to , and it was verified numerically that 
yield equal results. The integrals in Eqs. (4.26) and 
(4.28) may be considered as a higher order moment of the 
former integrals, hence, it becomes plausible to suppose 
that they should exhibit similar features. Therefore, one 
would be able to write 



-Ec(cyl) = 7r L Aoc = e e L , 



(4.29) 



and 



£ c (sph) = 4 7Ti? (sph) Aoc = a c N 1/3 . (4.30) 



Here Aoo is twice the asymptotic value of the integrals 
in Eqs. (4.26) and (4.28) and may be interpreted as the 
asymptotic coefficient for the curvature energy. It is wor- 
thy of notice that the curvature energy may be expressed 
in terms of the area and curvature of the cylindrical and 
spherical surfaces. According to a theorem of Euler the 
average curvature at a given point of a surface is defined 



as' 



32 



C- - (— — 
~ 2 \~Ri~ + ~R~ 2 



(4.31) 



where R\ and i? 2 are the radii of curvature along any two 
orthogonal tangents. Hence, for wide cylinders one gets 



1 



1 / J_ 

cyl ~ 2 \Kp + Rz -> oo; ~ 2i? (cyl) ' 



1 



for large spheres 



1 / 1 



-'sph 



2 \R 



1 

Re 



1 



i?o(sph) 



(4.32) 



(4.33) 



while for planar slabs it holds C s i aD = 0. So, one may 
write 



6 



E c (cy\) = 2TrR (cy\)L 



1 



2i?o(cyl) 



Aqc — ^cyl ^cyl Ao 



(4.34) 



and 

E c (sph) 



4 7 ri? 2 1 (sph) 



-Ro(sph) 



Xr 



^sph ^sph 



(4.35) 



The relation between a c and et may be derived from the 
ratio of these equations 



-E c (sph) _ AjphCsphAoc _ 4i? (sph) 



^c(cyl) 
which yields 

leading to 



^4cyl ^cyl -^oc 



4fl (sph) 

JVl/3 



El L 



47T po 



1/3 



48 / e £ ~ ' 



(4.36) 



(4.37) 



(4.38) 



The values of this relation calculated with parameters 
determined from both examined DF arc included in Table 
I. For the Skyrme-DF the result was 0.99 in excellent 
agreement with the DM prediction. The ratio yielded 
by the OP-NLDF, i.e. 1.10, differs slightly from unity. 
This is indeed to be expected when the DF is not only 
written in terms of p(r) but has explicit information on r 
as happens with Lennard-Joncs potential. In such a case 
there is an additional curvature contribution which in 
the approach adopted here is embedded in the coefficient 
a c - 31 

Upon starting from Eqs. (4.11) and (4.18), in both 
cases, one may express the excess of surface energy in 
terms of the curvature and the quantity Aoo getting 



(o/l 



\ Pa -i/2 

Coojcyl - \J — £ e n \ 



and 



(a a - Coo) S ph — a c 

'-)' 



7ri? (cyl) 

2 Aqc ^cyl j 
1/3 



(4.39) 



2pjT) at-i/3 
9tt 



/it p \ V 3 
V 6 



1/3 



' ^sph 



7T i?o( s P n ) 

= 2 Aoo C sph . (4.40) 
By taking into account these results one arrives at 

{a A - Coo)sph Asph _ ^sphCsph _ -E c (sph) 



oo Jcyl ^-cyl 



j4 cv i C, 



cyl ^cyl 



E c {cy\) 



(4.41) 



This means that the ratio of the excess of surface energies 
is equal to the ratio of curvature energies. 



B. Thickness at the free interface 

All the calculated results for the surface thickness W 
of free 4 He systems are plotted together in Fig. 5 as a 
function of e — es- Let us first look at values obtained 
by applying the Skyrme-DF. In this case we must state 
a word of caution because our results for spheres differ 
noticeably from those listed in the fifth column of Tabic 
I in Rcf. 13. We attribute this difference to a misprint in 
Table I of Rcf. 13 in view of the fact that the tabulated 
values are also inconsistent with the thickness inferred 
from the profiles plotted in Fig. 1 of the same paper. 
This point is stressed because due to that mistake the 
range of W ascribed to calculations with the Skyrme-DF 
for spheres with N < 728 quoted in the summary of Table 
III in Rcf. 4 is incorrect, actually itis6.6<I4 / <7.4A. 

Turning to the present results, the displayed data indi- 
cate that the asymptotic surface width obtained with the 
zero-range functional does not depend on the geometry of 
the helium system. In all the cases it attains Woo - 7 A, 
in agreement with the asymptotic result for an infinite 
system reported in Rcf. 13. 

The results obtained from calculations carried out with 
the OP-NLDF are smaller than those yielded by the 
Skyrme-DF, but exhibit a similar trend as a function of 
e— es and the asymptotic width is also independent of the 
geometry being ~ 5.9 A. For large systems the OP- 
NLDF thickness is close to that reported for droplets in 
Table II of Dalfovo et al. 4 These authors have performed 
calculations for large systems (1000 < N < 10000, i.e., 
-6.4 < e < -5.4 K) with the more elaborated OT-NLDF 
getting W ~ 5.7 A, as may be seen in the drawing. 

The values obtained from Monte Carlo simulations for 
free planar slabs 6 ' 12 and spherical drops 6 are also plotted 
in Fig. 5. However, in particular, the data evaluated 
by using glue-SWF 6 exhibit a too large spread to make 
possible any meaningful systematics. 

The comparison between the theoretical results men- 
tioned above and the experimental data is performed in 
the same Fig. 5. The data obtained in the measurement 
of large droplets published in Ref. 4 are represented by 
the reported mean thickness W = 6.4± 1.3 A. In the case 
of Ref. 7 only the two data included by the authors in 
the abstract are shown. Since these values correspond to 
rather broad planar films, we plotted them schematically 
close to origin for the abscissa. A similar criterion was 
adopted for the value of Ref. 3. As can be seen in this 
figure, all the results provided by finite range function- 
als for the largest systems are consistent with the most 
recent and precise experimental data. 7 

The measurements of Refs. 3 and 7 have been per- 
formed for planar films of 4 He adsorbed onto a substrate 
of Si, while the calculations mentioned above have been 
carried out for free planar slabs. Therefore we completed 
the investigation exploring how much theoretical results 
change when helium is adsorbed on a solid surface. In 
doing so, we include in this analysis the results of our 
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OP-NLDF estimations of the thickness at the free in- 
terface for 4 He films adsorbed on surfaces of Na and Li. 
These substrates provide the strongest adsorption poten- 
tials among alkali metals. 33 ' 34 For example, we show in 
Fig. 5 the values calculated for stable films with cover- 
age n c = 0.34 A~ 2 in the case of both substrates, i.e., 
W = 5.57 (Na) and 5.45 A (Li). 35 These results fit per- 
fectly into the general pattern exhibited by evaluations 
performed with NLDF and recent experimental data. 7 



1. Simple models for the surface thickness 

Let us now show that very simple models provide rea- 
sonable estimations of the interface thickness and its re- 
lation with the surface tension for large systems. For this 
purpose, it is useful to start first from the Eq. (2.6) of Ref. 
5 for a a of planar slabs obtained within the Skyrmc-DF. 
It reads 



OA 




2d 4 



dp{z) 
dz 



(4.42) 



We shall evaluate this integral adopting a crude approxi- 
mation for the density profile of a wide slab. Let us con- 
sider that from z = to z = zq the density is constant 
and equal to its value at the center p c . Subsequently, 
for z > z begins the falloff. Assuming that the first 
stretch of the falloff up to z = Z£ is linear, provided that 
p(zt) < /3 C /10, one may write 



p{z) = Pc 



1 - 



w 



4 (Z - Zq) 

5W 



(4.43) 



Next, for z > zg, the falloff of the density profile follows 
the asymptotic exponential law (see, for instance, Eqs. 
(7) and (9) in Ref. 36) 



p{z) 



a 



exp [—2 B (z — ze)] 



(4.44) 



with B = y — 2mes/?i 2 = 1.09 A 1 . The parameters zi 
and a are determined by matching at z = zt the expres- 
sions (4.43) and (4.44) for p(z) as well as its derivatives 
dp(z)/dz. This task is performed in Appendix B and the 
result is 



OA 



Pc 

5 W 



+ 8d 4 p c 



(4.45) 



So, in this crude approach the asymptotic width be- 
comes 



Woo ^^(-+8d 4 p 
5 (Too V rn 



6.83 A . 



(4.46) 



lying close to the limit suggested by data displayed in 
Fig. 5. 

A slightly different version of the simple model de- 
scribed above may provide separate formulas for the sur- 
face thickness and tension in terms of parameters of the 
Skyrmc-DF. As we shall see below, such expressions are 
valid for all the analyzed geometries. In this model, 
the proposed p(r) is the simplest that take into account 
three of the most relevant characteristics of density pro- 
file, the sharp radius of the density distribution R s , the 
thickness of the interface W, and the "leptodermous" 
behavior p c > po- We assume that from r = to 
r = Ri = R s — (3W the density is constant and equal 
to p c = p(r = 0), we recall that W = 5W/4. At r = Ri 
begins a linear falloff and for r > i?2 = Rs + (1 — P)W 
the density vanishes. The parameter (3 depends on the 
geometry through the normalization condition given by 
Eq. (2.3), where N should be suitably written for non 
spherical systems in terms of n c or n\. The obtained for- 
mulas for the normalization condition and the parameter 
(3 are given in Appendix C. Subsequently, the ground 
state energy provided by the Skyrme-DF for slabs, cylin- 
ders, and spheres was expressed by taking into account 
the proposed p(r). In order to avoid an unphysical di- 
vergence in the integration of the kinetic term caused by 
the kink of the adopted density profile we assumed 



dz 



5 W 



(4.47) 



As an example, the energy for cylinders is given in Ap- 
pendix C. 

Upon minimizing the grand free energy (e.g., for cylin- 
ders see Eq. (C5)) with respect to W and e, and keeping 
only the dominant terms, the following relation valid for 
all three considered geometries is obtained 



Wa A = ^(-+2d 4Pc 
5 V m 



(4.48) 



Here, the dominant contribution which is given by the 
term proportional to d 4 is equal to that previously ob- 
tained in Eq. (4.45). This procedure also yields an ex- 
pression for the asymptotic width in terms of parameters 
of the functional 



VFoo = 



yj (h 2 /m) +2d 4 p 

PB T 3 PO "I- 74+3 Po 



= 6. 90 A. (4.49) 



This value is in excellent agreement with the asymptotic 
result obtained from the complete numerical solutions in 
the Skyrmc-DF approach. The corresponding result for 
the asymptotic surface tension 

Oqo ~ TtT" - ( — + 2 d 4Po ] = 0.294 K/A 2 , (4.50) 
« " oo \rn J 

indicates a good degree of self-consistency of this approx- 
imation. 
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In conclusion, we can state that the crude simple mod- 
els used here to estimate the surface thickness are able 
to catch a substantial part of the physics involved. 



2. The quantity Wry a 

By looking at the relation between W and a a given 
by Eqs. (4.45) and (4.48), we found it worthy to examine 
the product W a a as a function of e — eg. Such a behav- 
ior is shown in Fig. 6, where for curved geometries one 
may observe overshoots similar to that exhibited by the 
central density and the surface tension. Since according 
to Eqs. (4.45) and (4.48) the main contribution to W a a 
is given by the term containing di, in this case the linear 
departure would be governed by 



p 2 c =p 2 (l + e) 2 ~p 2 (l + 2e), 



which leads to 



Wa A ^(Wa A )oo(l + 2e) 



(4.51) 



(4.52) 



This relation yields linear asymptotic expressions, where 
the slope is determined by the compressibility /C which 
is an input for the adopted DF approaches: 
a) for cylindrical systems 



W<t a ^ {Wa) z 



b) for spherical systems 



W a a — {Wa) c 



(4.53) 



(4.54) 



These approximations are plotted as dashed lines in Fig. 
6. From this drawing one can realize that these very sim- 
ple forms account fairly well for the obtained departures. 
In the case of planar systems there is no overshoot. 

By combining the expansions for the product Wua 
and for the surface tension a a it is possible to estimate 
the departure of W from its asymptotic value: 

a) for cylindrical systems 



OA 



1 + 



Po £e 



(e - e B ) 



(4.55) 



b) for spherical systems 
{W <ta) 



W 



OA 



4 



Pa a c 



n Pa 
6 



1/3 



(4.56) 



For both geometries there is an important cancellation 
within the coefficient of the linear term in these expres- 
sions. This fact explains why for large systems the thick- 
ness as a function of e — is flat in Fig. 5. 



V. SUMMARY 

The curvature effects on the thickness of the free in- 
terface W and the corresponding surface tension a a are 
studied. It is known that the free spherical systems are 
those with lowest surface area to volume ratio which leads 
to lowest energy per particle, but nevertheless, we found 
of interest to examine also other geometries to have a 
more complete view of surface properties. In addition, 
some general properties of free systems could be extended 
to adsorbed ones. 

First of all, we tested whether the largest systems have 
already reached the critical size for exhibiting a global 
asymptotic behavior. From Figs. 1 and 2 one may real- 
ize that the moments < r k >, with k — 1 and 2, have 
attained the asymptotic values. 

In this paper, the results for a a and W corresponding 
to different geometries are presented in a unified scale as 
a function of e — es allowing a direct comparison. As far 
as we know, this procedure has not been previously used. 

It was found that for cylindrical and spherical systems 
the surface tension presents an overshoot similar to that 
previously observed for the central density. 13 ' 18 This fact 
indicates that the squeezing effect known as the "lepto- 
dermous" behavior is also manifested in the surface ten- 
sion. On the contrary, this feature is not present in the 
case of symmetric planar slabs. It is shown that the be- 
havior of large systems may be well interpreted within 
the DM 30,31 adapted to each geometry. In particular, a 
relation between coefficients of the curvature energy is 
derived. 

An analysis of the surface thickness W indicates that 
for large systems, i.e., when e — < 1 K, the asymptotic 
result is independent of the geometry of the 4 He system. 
The departure from that value is different for each ge- 
ometry. These features are common to both examined 
density functionals. 

We should mention that according to a discussion given 
in the text, the range 8.8<VF<9.2A quoted in Table 
III of Ref. 4 as corresponding to droplets (with N < 728) 
studied by utilizing the Skyrmc-DF, which was deter- 
mined with data taken from Table I of Ref. 13, must be 
replaced by 6.6 < W < 7.4 A. 

The results obtained with the OP-NLDF for large sys- 
tems are close to those calculated with the OT-NLDF for 
spheres by Dalfovo et al. 4 It is shown in Fig. 5 that the 
most recent experimental data 7 favor the results yielded 
by the finite range density functionals. In addition, the 
thickness of the free surface of helium films adsorbed on 
planar substrates of Na and Li evaluated with the OP- 
NLDF matches very well with the pattern depicted in 
Fig. 5. 

We have also analyzed the product Wga- As expected, 
for curved geometries this quantity exhibit an overshoot. 
However, it is interesting to note that in this case the de- 
parture from the asymptotic values is mainly determined 
by the inverse of the compressibility JC of bulk 4 He. Fur- 
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thermore, this analysis gives a hint for understanding the 
flat behavior of W at small e — e#. 
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APPENDIX: A. OP-NLDF HARTREE 
MEAN-FIELD POTENTIAL 

In this Appendix we compile the relevant expressions 
derived in the OP-NLDF approach for a spherical geome- 
try. The explicit form of the Hartree mean-field potential 
Vsr(r) derived according to Eq. (2.5) becomes 

+y(74 + l) J drfp(r/)[p(rf)Y l4 W(\ r-r/ |) . (Al) 

Let us first provide the expressions of the contributions 
involving the "coarse-grained density" p(r), i.e., 

p{r) = j drl p{ri) W{ | r - r/ |) , (A2) 

and 

pv{r) = j dr/p(r/)[j5(r/)}^W(\r-r/\) . (A3) 
Both these integrals may be cast into the form 



H(r) 



Airh 3 0¥ 



dr/TZ(r/)Q(h p- \ r-rl |) . (A4) 



After introducing spherical coordinates and taking into 
account that the step function is symmetric in the az- 
imuthal angle ip, the integration over this variable yields 



j rl 2 drlK{rl) j 
x 0[h 2 OP - r 2 - r/ 2 + 2rrf cos0] . (A5) 



K{r) = — — I rl' 2 drm{ri) \ sin 0d6 

2 O.QP Jrl, 



For the point located at r = one gets 

O /''''max 

n ( r = °) = oTs~ / r ' 2 drl U ^ 6 1 ^op - r ' 2 ] 
2 h OP J r/min 

x r 

sin 9 d6 



r ■ 

/ sir 

Jo 

I 



/iqp 



r/ 2 dr/TZ(r/) . 



(A6) 



For r > two different cases should be considered: (i) 
< r < h G p and (ii) r > h Q p. For < r < h p the 
integral over rl should be split into two parts 



H(r) = 



3 



2 ft-QP i- Jo 



h p-r 



rr+hop ~j 
Jhov—r ) 



ri 2 drf1Z(rf) 



f 

Jo 



sin d9 9 [ h 2 ov - r 2 - rl 2 + 2rr/ cos ] . ( A7) 



Since for the first integral over rl the upper angular limit 
is #max = 7T 5 while for the second integral # max is deter- 
mined by the condition 



cos n 



r 2 + rl 2 - h 2 



OP 



2rrl 



(A8) 



then 
H{r) 



9 h 3 

z ' 1 op Jo 



hop-r 



ri 2 driTZ(ri) / sin0d0 



+ 



r+h F 



Oh 3 

^ n OP Jhop-r 



rl 2 drl 



1l{ri) f 
Jo 



sin9d0. (A9) 



The integration over leads to a very simple expression 



H(r) 



h 3 

,l op Jo 



h p-r 



ri 2 dr/lZ(ri) 



+ 



J 

Jh 



P Jhop- 



ri drfK(rf) 



r — rl 
hop 



(A10) 



For r > hop, there is only one contribution similar to 
the second term in Eq. (A9) 



TZ(r) 



9 /i 3 

z ' l OP Jr-hov 



r+hop 



ri 2 driK(ri) 



4r/i OP J r - ho 



r+h p 



r/dr/1l(r/) 



sin0d0 



r — rl 



hoi 



(All) 



Let us now focus on the integration of the screened LJ 
potential contributing to Eq. (Al) 



V-LJScr 
H 



(r) 



dr/p(ri) Vj OP (| r rl |)) 



Jr. 



with 



+ 4e / dr/p(ri) 

lR>h p 



R=\r-ri 



(iT-(i) 



, (A12) 



(A13) 



Here a and e stand for ctl.t and clj, respectively. In this 
case one can follow the same procedure as that utilized to 
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calculate the "coarse-grained density" terms. As before 
there are three different domains of r to be considered. 
At r = the Hartree potential reads 



F LJScr( r = Q ) 



AlT 



r^OF 



+ 16tt£ / rl 2 drl p{rl) 
J hop 



f 

M Vr P (h OP ) ri 6 drfp(ri) 
n OP Jo 



(- 

\ri 



(A14) 



For < r < hop the integral over r/ becomes 



V^ ,bcI {r) 



'OP TfOP 



^ uf (^op) 



3r 

hop—r 

x< / rldrlp{rl) 



r + rl 
hop 



r — rl 
hop 



6n 



+ ■ 



47T6 (7 



rr+hop 

+ / rl drl p(rl) 

Jhop—r 

2 r r + h op 

rl drl p(ri) 

hop—r 
10 



+ ■ 



„ r + r/ _ 



rl drl p(ri) 



r+hop 
10n 



1 - 
1 

5 



r + rl 



r — rl 



r — rl 
h p 
a 
hop 
a 

r + rl 
a 

r — rl 
a 



10 



10 



r + rl 



(A15) 



Finally, for r > hop one gets 
V L H 1S "(r) = ^V l OP (hop) 



3r 



rr+hop 

/ P{r') 

Jr—hap 



r — rl 



'OP 



rl drl 



| Aire a 2 f r' ,OP + 



r — rl 
a 

r — rl 



10 



r/ dr/ p(r/) 

10n 



r + r/ 

<7 



r + rl 



Aire a 2 f r+hop , , , f 1 
H / rl drl p{rl) 

^ Jr—hop 



( a ) 


1 


\r + r/J 


~ 2 



ft-OP 



o 
hop 
a 



r + rl 



(A16) 



P(^) = Pc ( 1 - 

and 

dpjz) 
dz / 

These conditions lead to 



4(££--£o) 

51f 



4p e 
5VF 



a 



a = - BVF 



and 

ze - z o 



■W 



1 
2~B 



W- 



w 



1 

~2B 



(Bl) 



(B2) 



(B3) 



(B4) 



Since for large slabs one expects W > A/B ~ 3.7 A, 
then zi — z would be larger than 9 W/8 ensuring p(zg) < 
p c /10, which in turn would support the consistency of 
the adopted law for the falloff and the definition of the 
thickness. 

For the evaluation of the derivative one may use the 
relation 



d^fpjz) 



1 



dz 



dp{z) 

2^/pJz) dz 



(B5) 



Now, the integral of Eq. (4.42) may be split into two 
parts 



a a = 



which leads to 

Pc 



dz 



1 



h z 



Am p{z) 



2d 4 



dp{z) 
dz 



(B6) 



° A 5 W I m 



1 + In ( - B W 



■8 C?4 p c 



1 



1 



5BW 



(B7) 



It is possible to verify that for the expected values of W 
there is an important cancellation between terms carry- 
ing the product B W. Hence, Eq. (B7) may be reduced 
to 



° A - tttt — v 8 di P c 



APPENDIX: C. SIMPLE MODEL II 



(B8) 



APPENDIX: B. SIMPLE MODEL I 

Starting from Eqs. (4.43) and (4.44), after imposing 
continuity of functions and derivatives one gets 



The evaluation of the normalization condition yields 
[R, + (1 - [3)W} D+1 - [R s - (3W} D+1 



n = n D Pc 



D(D + 1)L D - 3 W 



(CI) 
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with Qd=i = 2, il_D = 2 = 2 7r, and £7 = 4 7r. In the 
case D = 1 the "sharp radius" i? s stands for z s . The 
parameter /? becomes: 

a) for the planar geometry 



= 



1 

2 5 



b) for cylinders 

c) for spheres 



' 2 W 2V3 





+ - 



2 l/3 

12~ 



+ 



2 * 6 3 



-1/3 



(C2) 



(C3) 



1/3 



(C4) 



For instance, the ground-state energy for cylindric ge- 
ometry reads 

e = e B + 21m- 1 1 l eB (3W(-2R s + f3W) Pc 

1 R s \ ( n 2 , \ 

2 + ^-^j U^H* 



3c4 



+■ 



1 



74 + 3 V74 + 4 W 



(C5) 



Here p c and i? s are functions of e given by Eqs. (4.1) and 
(4.7), respectively. 
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TABLE I. Bulk observables for liquid 4 He at T = and 
the calculated parameters ee and a c . 



Observable 


Data 




Ref. 


e B [K] 
Po [A- 3 ] 
K [K] 

a cxp [K/A 2 ] 


-7.15 
0.021836 
27.2 
0.274 ± 0.003 
0.257 ±0.001 
0.272 ± 0.002 




21 
21 
21 
8 
9 
10 


Parameter 


Value 


Theory 


Ref. 


e e [K/A] 
£e ["/ aj 
a c [K] 
a c [K] 
a c [K] 
a c [K] 


1.237 

U.OOi 

10.45 
10.90 
10.86 
8.58 


Skyrme-DF 
DP NT DP 
Skyrme-DF 

OP-NLDF 


PW a 

PW 

L VV 

13 
21 
PW 
PW 


( n p \ V 3 a c 

V 48 J e e 
I it p a \V 3 a c 

V 48 J e t 


0.99 
1.10 


Skyrme-DF 
OP-NLDF 


PW 
PW 


( it po \ V 3 a e 
V 48 1 e e 


1.00 


DM 


Eq.(4.38)-PW 



a PW stands for results obtained in the present work. 
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FIG. 1. Normalized moments [< r k (n x ) >] 1/fc /n^ /2 as a 
function of the inverse of square root of longitudinal density 
v cy \ = n ^ 1//2 f° r f ree helium cylinders. Open and full sym- 
bols stand for Skyrme-DF and OP-NLDF results, respectively. 
The dashed and dot-dashed curves are only given to guide the 
eye. The solid lines correspond to the asymptotic values given 
byEq. (3.5). 

FIG. 2. Similar to Fig. 1, but for free helium spheres. The 
solid lines correspond to the asymptotic values given by Eq. 
(3.7). 

FIG. 3. Surface tension calculated by using the Skyrme-DF 
is shown as a function of the energy difference e — es. In 
order to not overcrowd the picture only some selected data are 
plotted and solid curves are drawn to guide the eye. Dashed 
lines indicate the linear asymptotic laws for curved geometries 
given by Eqs. (4.21) and (4.22). 

FIG. 4. Same as Fig. 3 with data calculated by using the 
OP-NLDF. In addition, the dashed curve for planar slabs is 
the quadratic asymptotic law given by Eq. (4.23). 

FIG. 5. Thickness of the free surface as a function of the 
energy difference e — en- The circles, triangles, and squares 
represent own data for spherical, cylindrical, and planar sys- 
tems, respectively. Empty symbols stand for data calculated 
by using the Skyrme-DF, while full symbols are results from 
the OP-NLDF approach. As in previous figures, only some se- 
lected data are plotted and solid curves are drawn to guide the 
eye. Diamonds are OT-NLDF results for spheres from Ref. 4. 
Full and empty five points stars are Monte Carlo simulations 
for spheres and slabs reported in Ref. 6. Six points stars are 
Monte Carlo data from Ref. 12. Asterisks are OP-NLDF re- 
sults for helium adsorbed onto planar Na and Li substrates. 
The cross (x), the "otimes" ((g)), and the "oplus" (©) are 
experimental values from Refs. 3, 4, and 7, respectively. 

FIG. 6. Quantity Wo a as a function of the energy differ- 
ence e — es. As in Fig. 5, the circles, triangles, and squares 
represent data for spherical, cylindrical, and planar systems, 
respectively. Empty symbols stand for data calculated by us- 
ing the Skyrme-DF, while full symbols are results from the 
OP-NLDF approach. Dashed lines indicate the linear asymp- 
totic laws for curved geometries given by Eqs. (4.53) and 
(4.54). 
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